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Abstract. We study discretisation effects in cellular automata models for pedestrian 
dynamics by reducing the cell size. Then a particle occupies more than one cell which 
leads to subtle effects in the dynamics, e.g. non-local conflict situations. Results 
from computer simulations of the floor field model are compared with empirical 
findings. Furthermore the influence of increasing the maximal walking speed v max is 
investigated by increasing the interaction range beyond nearest neighbour interactions. 
The extension of the model to w max > 1 turns out to be a severe challenge which can 
be solved in different ways. Four major variants are discussed that take into account 
different dynamical aspects. The variation of ^ ma x has strong influence on the shape of 
the flow-density relation. We show that walking speeds v max > 1 lead to results which 
are in very good agreement with empirical data. 
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1. Introduction 



The understanding of the dynamical features of pedestrian dynamics has been the aim 
of many investigations over the last years PjJ 12 E|. Recently also approaches based 
on cellular automata (CA) have been suggested. CA models are discrete in space, 
time and state variables. On the one hand, this makes the models ideally suited 
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for large-scale computer simulations. On the other hand, the discreteness has to be 
regarded as an approximation of reality. In this paper we want to investigate some more 
fundamental questions related to general aspects of physical modelling using discrete 
models. More specifically we try to elucidate the qualitative and, in some cases, 
even quantitative influence of the parameters related to the discretisation. As with 
any other kind of simulation (e.g. fluid dynamics j3J), two major decisions have to be 
made for the description to be used: macro- vs. microscopic and discrete vs. continuous 
[Sj. Although we do not consider it explicitly here, performing the continuum limit 
of the discrete models would contribute to a better understanding of their relation to 
continuous models. 

As we will see, the discreteness leads to some problems that will be exemplified 
using the so-called floor field CA model that has been introduced in (fHHHE]- The big 
advantage of this model compared to other CA approaches for pedestrian dynamics 
j3 EDI HU 112 EE] is that the floor field model is able to reproduce most of the 
characteristic aspects of pedestrian dynamics, especially the different collective effects 
observed empirically PjJ I2EJ E3- Here, the inclusion of friction effects allows to describe 
many important observations, especially in egress and evacuation scenarios [HI Elj • 

There are foremost three basic parameters any CA model is based on: the 
interaction range, which in traffic models usually is given by the maximal velocity 
v max , the time scale At and the generic length scale a, corresponding to the size of 
one grid cell. The variation of At will only cause the rescaling of all time values such 
as averaged evacuation times and is therefore from a theoretical point of view not very 
interesting. In contrast, variations of t> max and a will turn out to have a strong qualitative 
influence on the dynamical properties of the model. 

The outline of the paper is as follows: First, we present extensions of the basic 
model (with t> max = 1) to larger interaction ranges t> max > 1. Four possible variants of 
the extension are considered, based on different aspects of the dynamics, especially the 
treatment of crossing paths. Monte Carlo simulations will be presented and compared 
to empirical data. Then we investigate the effects of a reduction of the length scale a. 
We will also show that a finer spatial discretisation (by reducing a from a = 40 cm to 
a = 20 cm) has strong effects on the dynamics. 

But first, we will review the most important characteristics of the basic model with 
Vx = 1 and a = 40 cm in the next section. 

2. Cellular automata models 

For pedestrian dynamics, in a natural spatial discretization one introduces cells of 
40 cm x 40 cm. This corresponds to the typical space occupied by a person in a dense 
crowd. Therefore, each cell can then either be empty or occupied by (at most) one 
pedestrian (hard-core exclusion). 

Most CA approaches use a parallel (synchronous) update where the dynamical rules 
Z have to be applied to all particles at the same time. This introduces a timescale At. 
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In stochastic CA models the motion to other cells at each discrete time step t — > t + At 
is controlled by transition probabilities. In the simplest case, which is used most often, 
each pedestrian is allowed to move only to one of his/her four nearest neighbour cells. 

The effects of the discretisation and interaction range that we want to investigate 
here will appear in any CA model. However, our testing ground will be a specific CA, 
namely the floor field model. It can be regarded as a two-dimensional generalization of 
the asymmetric simple exclusion process (ASEP) (see e.g. where the transition 

probabilities are determined dynamically through a coupling to so-called floor fields |EllZj- 
In turn, these floor fields are changed by the motion of the pedestrians which leads to 
feedback effects. This way of implementing the interactions is inspired by the form of 
communication in insects societies (e.g. ants) and the interaction can be described as 
virtual chemotaxis. Generically a movement in the direction of larger fields is preferred. 
The key point is that an interaction that is long ranged in space can be translated into 
a local interaction, but with "memory" . Note that in contrast to chemotaxis and trail 
formation [2131211221123] in our model the trace is only virtual, i.e. can not be observed 
empirically. We do not give a full definition and discussion of the floor field model, 
which can be found in [HI El 123 EH] but focus on the aspects that are important in 
the following. 

The floor field consists of two parts, a static and a dynamic field. Both are discrete 
and we can imagine the field strength as the number of virtual pheromones present at a 
site. The fixed static floor field S is not changed by the presence of the pedestrians. S is 
used to specify regions of space which are more attractive, e.g. exits or shop windows. 
In case of an evacuation simulation the static floor field S describes the shortest distance 
to the closest exit door measured in steps, i.e. the number of lattice sites to the exit 
[7j. Steps can be carried out to the four next neighbour sites. The field values increase 
towards the door. 

Collective effects like lane formation or herding PJ 121 ESI EZJ that are based on a 
long range interaction between pedestrians can be taken into account by a dynamically 
varying floor field D H- The dynamic floor field D is a virtual trace left by the 
pedestrians and has its own dynamics, i.e. diffusion and decay. It is used to model 
different forms of interaction between the pedestrians. At t = the dynamic field is set 
to zero for all sites (x,y) of the lattice, i.e. D xy = 0. Whenever a particle jumps from 
site (x, y) to one of the four neighbouring cells, D at the origin cell is increased by one: 
D xy — > D xy + 1. Furthermore, the dynamic floor field is time dependent, it has diffusion 
and decay controlled by two parameters < a, 5 < 1, which leads to broadening and 
dilution of the trace [Hj . 

The coupling of the particles to the static and dynamic floor field is controlled by 
coupling constants k$ and k D , respectively. A strong coupling to the static field implies 
an almost deterministic motion in the direction of larger fields. On the other hand, if 
the coupling to the dynamic field dominates, herding effects become important where 
pedestrians blindly follow others as it may happen in the case of panic. 

An important extension that makes the dynamics more realistic is the concept of 
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friction [SJ EI] • It has been implemented in other models in a direct way as contact 
friction [27J. However, in our approach it concerns the behavior of two (or more) 
pedestrians that try to move to the same cell. With some probability /i, called friction 
parameter, none of the particles is allowed to move, whereas with probability 1 — \x one 
of the particles is chosen randomly to occupy the target cell. This leads to effects similar 
to arching in granular materials [27[ I2E] • 

3. Extensions to f max > 1 

In the definition of the model given above, particles are only allowed to move to 
unoccupied nearest neighbour sites in one time step t — > t + At, i.e. t> max = 1 cell/time 
step. This constraint corresponds to a maximum empirical walking speed of v emp = 1.3 
m/s [HI HHj, if the time scale At is identified with 0.3 s. The choice t> max = 
1 has the advantage of simplicity and high simulation speed. It also reproduces 
most experimentally observable phenomena. Nevertheless, there are reasons for the 
introduction of higher walking speeds t> max > 1 in some cases: 

(i) Fig. E shows the fundamental diagram of the flow along a corridor in the v max = 1 
case of large coupling k$ and hp = 0. § It is nearly symmetric with maximal 
flow at p ~ However, experimental data point to a non-symmetric fundamental 
diagram with maximal values of the flow for densities p < \ E3 • 

(ii) The interaction horizon is not isotropic since pedestrians react mainly to stimuli in 
front of them. This anisotropy can be better taken into account in a model with 
larger interaction range where it is generated dynamically through the exclusion 
principle. 

(iii) Although the velocity distribution of a pedestrian is sharply peaked around 1.3 m/ s, 
higher walking speeds are frequently observed [TU HH] . 

(iv) In a model with f max > 1, a realistic distribution of different walking speeds of the 
particle ensemble can be implemented in a simple way. 

Certainly the most relevant aspects are the first two. The other two requirements could 
be taken into account in the case f max = 1 by introducing individual randomization 
probabilities, determined by the coupling constants ks and ko, or a rescaling of the 
relevant time scale. In the following extensions of the model to t> max > 1 will be 
presented. The possible variants of the update algorithm and their dynamical properties 
will be discussed in detail and compared to experimental data. 

3.1. Classification of the model variants 

In the case of f max = 1, the particles are allowed to cover a maximal distance of a = 40 
cm (corresponding to a movement to next neighbouring sites) per time step. Then 

§ The empirical unit of the specific flow is pedestrians/(mctcr-sccond), the natural unit of the model: 
particles/ (a • At) (with cell length a =40 cm and time scale At). 
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Figure 1. Fundamental diagram (flow J, measured in units of 1/At, vs. density p, 
measured in units of p max = l/a 2 = 6.25 m~ 2 )) for a corridor of 93 x 33 lattice cells 
of size a 2 — 40 x 40 cm 2 with i> ma x = 1 in the so-called ordered regime (ks = 10.0 
and kn — 0.0). The static floor field S is calculated using a Manhattan metric. In the 
direction of the flow periodic boundary conditions are used. For details see |7II21). 



for higher walking speeds w max > 2 the particles should be able to cover a distance of 
f; max ■ a in one time step. This does not necessarily correspond to higher velocities, e.g. 
when combined with a reduction of the cell size or change of time scale, which is the 
case mainly considered here. Then, due to the two-dimensional nature of the motion, 
the choosen path is not necessarily straight. Also it has to be taken into account that 
now the paths of different pedestrians may cross. This requires a generalization of the 
concept of friction. In the following we develop a formalism which allows to describe 
the different possible generalizations. 

At the beginning of a time step t the position of each particle n e {1, . . . , iV} is 
denoted by TQ U \t). Each particle then chooses, corresponding to its desired movement, a 
trajectory T^ n \t) = {TQ U \t), . . . , Tv2L(t)} consisting of v max + l lattice sites. The single 
elements of T^ n \t) are determined by applying the update rules Z for the t> max = 1 case 
successively t> max times. Here, every arbitrary element (t) (I > 0) corresponds to a 
target site Z of the particle n that has advanced virtually to site T^\(t) (see Fig. EJ: 



Therefore, 7] (t) can only be one of the four next neighbour sites of T^[(t) or 
T/ n) (t) =T/"j(t). Thus, every particle accesses iteratively at most v max target sites. It 
is important to notice that the v max target sites and thus the trajectories are determined 
using the configuration of the system at time t through {Tq (t), . . . , T (t) } . 



Tl n \t) = z(rn(t)). 



(1) 
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Figure 2. Graphical representation of the desired trajectory Tj (t) of pedestrian n 
for the case v max > 4. 



All desired trajectories can be summarized in the matrix 
( T (1) (t) ••• Til(t) \ 

T(t) = 



(2) 



which has to be stored at every time step of the simulation. Two basic assumptions of 
the model are reflected by the matrix elements of T(t): 

• Hard-core exclusion of particles: T^\t) ^ T^\t) for all i ^ j, i.e. every lattice site 
can be occupied by at most one particle. 

• Parallel updating: T^O) ^ T (i) (t) for all i^j. This expresses the reaction time 
of the particles. No particle can occupy a lattice site at time step t + 1, which was 
occupied at time step t by another one. 

These properties are special cases of the more general rule Tn\t) ^ T$\t) for all i ^ j 
and n — 0,1, ... , v max . They restrict the allowed values of the first and last column of 
T(t), i.e. the allowed positions at time t and t + At. The remaining t> max — 1 columns 
correspond to intermediate states that are not explicitly realized as configurations. 
It might happen that two or more desired trajectories cross or that even one of the 
intermediate states violates the hard-core exclusion principle. In the following we will 
denote such situations as conflicts. In principle one could ignore the existence of conflicts 
since they correspond only to virtual intermediate states. However, in order to make the 
dynamics more realistic it is necessary to take them seriously. Then a method for the 
resolution of conflicts is required. In the following we will discuss four major variants. 

3.1.1. Hop or stop The simplest variant ignores the intermediate positions T^ n \t) 
(1 = 1,... ,v max — 1) and takes only conflicts in the desired new positions T^^(t) into 
account. If k particles (k > 2) have chosen the same target cell, i.e. 

TSjt)=T^Jt), foralH,j'G{l,...,A;}, (3) 

one particle % is chosen by a probabilistic method (i.e. in random order) and is allowed 
to move to target site Ti^ ax (t) =: T^\t + 1). The other k — 1 particles j ^ i involved 
in the conflict are not allowed to move and have to remain at their origin sites, i.e. 
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T Uj (t+l) = T Uj (t). Practically, this can be easily implemented by allowing all particles 
to move sequentially by choosing randomly a permutation a of the particle indices j 
that determines the order a(j) of updating||. 

This is the simplest form of conflict resolution between particles. The overall 
movement takes place without consideration of the trajectories T^ n \t). Crossing 
trajectories and jumping over already advanced particles are permitted in order to keep 
the spirit of the parallel updating scheme. 

3.1.2. Move as far as possible This variant of conflict resolution is similar to the 
previous one. Again the particles are allowed to move in random order described by a 
randomly chosen permutation a. However, now other particles that have already reached 
their target sites cannot be jumped over by following particles anymore. A particle n is 
able to proceed along its trajectory T^(t), as long as T} (t) (I G {1, . . . , tWx}) is not 
already occupied by another particle. 

In this variant crossing trajectories are still possible, but jumps over particles that 
have already reached their target sites are not allowed. 

3.1.3. Vmax sub-steps In this variant of the update the particle movement along 
the trajectories T n (t) is subdivided into t> max sub-time steps where moves to nearest- 
neighbour cells are allowed. If in one sub-time step /, two or more (at most four) 
particles try to move to the same cell, i.e. 



only one randomly chosen particle is allowed to move forward to this site. All other 
particles then try to reach the site during the next sub-time step / + 1. 

In all sub-time steps of this update version conflicts between particles are taken into 
account and solved. Therefore, it is the easiest way to incorporate the concept of friction 
[Hj into the model with t> max > 1. Otherwise, this method leads to increased simulation 
times due to the successive treatment of the sub-time steps. Note, that the dynamics of 
this variant is different from the genuine f max = 1 -case, since the trajectories are 
fixed during the whole time step. 

3.1.4- No crossing paths In this last version, no particle is allowed to cross the 
trajectory of a particle that has already moved. The particles move in random order. 
If the path of particle n crosses that of a particle n' that has already moved (i.e. 
a(n r ) < a(n)), 



T, 



T, 



■Cj) 



for any i,j, 



(4) 



T, 



T^ n ■* , with I, k arbitrary, 



(5) 



then 




1 (1), ( T- 1 (2),..., ( T- 1 (7V). 



(6) 
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k is in this case assumed to be the smallest index in the trajectory of n' for which this 
condition holds. Therefore, in this variant a particle is allowed to move until it reaches 
a cell that is either occupied or has been passed by previous particle that has already 
been moved at the same time step. 

3.1.5. General aspects of the variants In all variants, at the end of each time step 
the field values of the dynamic field D at all lattice sites of the trajectories traversed 
by particles are increased by one. Afterwards D is modified by its decay and diffusion 
dynamics as described in [HI Ej . 

Each of the four update versions is stressing different aspects of the two-dimensional 
motion on the lattice. Fig. |3] shows a representation of all different methods for conflict 
resolution. All four variants reduce to the basic model described in Sec. El in the case of 
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Figure 3. Graphical representation of all variants of conflict resolution in the case 
Umax = 4. The trajectories T^ n \t) of the particles are shaded gray. The numbers of 
the particles correspond to the random order (given by the permutation a) in which 
the particles execute their steps. In variants 1), 2) and 4) this order is designated once. 
In the variant 3) one has to choose a new random order in every single sub-time step. 
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The characteristic time-scale At should not be influenced by the choice w max > 1- 
Therefore, for a given value of v max the corresponding maximum empirical walking speed 
can be calculated via 

= ^~=^-1.3m/ S (7) 

with a = 40 cm and At = 0.3 sec. 

In the following, the consequences of higher walking speeds on particle flows will 
be investigated and the differences between the four variants are discussed. Afterwards, 
the dependence of the evacuation time for a room on t> max is investigated. 



3.2. Simulations 

The particle flows J = p(v x ) presented in this section are measured in a corridor of 
size X x Y = 93 x 33 with periodic boundary conditions in x-direction. The motion 
in ^-direction is limited by walls. The static floor field S increases from left to right, 
i.e. S x ,y = x [TKj . We remind the reader, that in the floor field model a motion in 
the direction of increasing field strength is preferred. The lattice was initialized by 
randomly distributing N pedestrians corresponding to a density p = N/ (X • Y) . If not 
stated otherwise, the coupling hp to the dynamical field vanishes (k D = 0) and therefore 
the dynamical parameters of D can be chosen as a = and 5 = 1. 




Figure 4. Fundamental diagrams for v max = 1, ... ,5. a) Variant 4 (no crossing 
paths) with ks = 10.0 and kjj — 0.0; b) Variant 2 (move as far as possible) with 
k s = 2.0 and k D = 0.0. 

values J are shown for the case of a dominating coupling to the static field (ks 3> ko, 
corresponding to the ordered regime jjj). The update is based on the no crossing paths 
variant. The maximal flow is reached at smaller particle densities for growing f max . The 
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maximal flow is also increasing with v max . An unexpected property of this variant is that 
for densities p > 0.4 the flow J is largest in the case t> max = 1. This can be explained by 
the blockage of wide parts of the lattice, resulting from the non-traversable trajectories 
of the particles. The influence of this blockage increases with v max due to an increased 
'effective density'. The faster a pedestrian moves, the more cells she or he blocks. In this 
respect, the no crossing paths scheme leads to the highest effective densities (Fig. EJ). 

Because of the strong coupling to the static field S, which suppresses swaying in 
y-direction, the fundamental diagrams look very similar to those of the ASEP or other 
traffic models in the deterministic limit (Fig. If one chooses, like in Fig. for the 
second version of the update, a smaller coupling strength (ks = 2), the effective walking 
speed v e g of the particles is decreasing. For k s — > 00 the effective walking speed of a 
single pedestrian is v e g — > v mSuX , corresponding to a deterministic motion. On the other 
hand, for k s = the particle performs an isotropic random walk since the information 
about the preferred direction through the static floor field is not taken into accout. 
Therefore, for ks = 2.0 the shape of the fundamental diagrams is becoming smoother 
(Fig. 0Jd). The maximal values of J are smaller with a slight shift of the maximum to 
higher densities (for v mSuX > 1) and the slopes in the regions of small and high densities 
decrease. 

In the following, the consequences of the different conflict resolution schemes on the 
flow are discussed. Fig. shows the fundamental diagrams of all four variants for the 
velocities t> max = 2 and t> max = 4 and ks = 2.0. For t> max = 2 the differences between the 

a) b) 




Figure 5. Fundamental diagrams for all four update variants with ks = 2.0 and 
k D = 0.0: a) w max = 2; b) i> max = 4. 

variants are not strongly pronounced (Fig. EH)- But the flow values J are highest for 
the v max sub-steps and smallest for the no crossing paths scheme in all density regions. 

However, for t> max = 4 the differences between the variants are easy to spot. For 
particle densities p > 0.1 the no crossing paths scheme produces much smaller flow 
values than all other variants, because of the non-traversable trajectories which are four 
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sites long. The third version of the update (sub-steps) leads to the highest flows in all 
density regions since the resolution of conflicts in sub-time steps enables some particles 
to move forward in a later sub-time step. 

Finally, we compare the simulation results with empirical data. Fig. El shows a 
comparison of fundamental diagrams of the v max sub-steps variant (see Sec. I3.1.3j) with 
a fit of empirical flow-density relations [T3|. The data for the specific flow has been 
transformed to the units of the model. The fundamental diagram for t> max = 3 is in 




Figure 6. Comparison of fundamental diagrams of the third update version and the 
= 1 to u roax = 5 with experimental data from [14 (fc<j = 4.0 and 



velocities u n 



<D 



0.0). 



good agreement with the experimental data. In particular, the density values for the 
highest flow match each other very well. Since k$ = 4.0, the particles have a reduced 
effective walking speed v e g < f max - In the model, this results in a diversification of the 
flows for particle velocities between v max = 3 and lower values. For high density regions 
p > 0.6 the model yields flow values which are too small. The reason are lane changes of 
particles especially at higher densities. If the site in front of a particle is occupied, lane 
changes are possible. Such a changing does not increase the flow in positive x-direction, 
but often inhibits the movement of another particle (see Fig. EJ). The fundamental 



Figure 7. Decrease of the flow due to lane changes, 
diagrams of the other t> max values do not agree as well with the empirical data. 
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3. 3. Simulation of the evacuation of a room 

As a second scenario which helps to understand the effects of t> max > 1 we investigate 
its influence on evacuation times. The room to be evacuated is a square of 63 x 63 cells. 
The exit has a width of one cell and is located at the center of one wall. Initially the 
particles are distributed randomly and then try to leave the room due to the information 
they get through the floor fields. This scenario has been studied in detail in [7j where 
three different regimes ( ordered for ks ^> hp, disordered for fco» ks, and cooperative) 
have been identified. 

Fig. |H1 shows averaged evacuation times for different combinations of the coupling 
parameters ks and kp for the move as far as possible variant of the update. All other 
schemes show the same qualitative behaviour. The particle velocities are t> max = 1, . . . , 4. 
Fig. |Ht can be interpreted in the following way: If the coupling strength to the static 

a) b) 

I 1 1 1 1 5500 I 1 1 1 1 . 1 1 1 . 1 




Figure 8. Averaged evacuation times for a room of size 63 x 63 sites and an exit of one 
cell. The particle density at the beginning of the evacuation is p = 0.3 (1116 randomly 
distributed particles). Different i> max -values in the second version of the update (move 
as far as possible), a) kjj — 0.0; b) ks = 0.4 and 6 = 0.3 and a = 0.0 as parameters 
for the dynamics of D. 

field S is very small (ks <C 1), the particle motion resembles a random walk. But even 
if there is only slightly directed movement, the variance in lateral direction is increased 
with f max . Therefore, the particles can reach the exit by chance faster for increased v max , 
resulting in decreased evacuation times. In the limit ks — > oo the particles move on the 
shortest possible way to the exit, the movement becomes deterministic. For higher 
particle densities just after the evacuation starts a big jam forms in front of the exit. 
Therefore, the averaged evacuation times are in this case determined by the particle 
number (via the flow through the door) rather than the maximum walking speed. Thus, 
for higher densities and ks — > oo the evacuation times using the move as far as possible 
strategy (see Fig. |H]top) are the same for all t> max values. 

A coupling strength of ks ^> 1 is presumably most relevant for a comparison with 
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real pedestrian behaviour. Therefore, for many applications the restriction to t> max = 1 
is adequate to reproduce realistic behaviour. In other cases (see above) a choice of 
■^max > 1 is indispensable. 

Finally, Fig. [HJd shows the different evacuation times for kp = 0.0 and varying f max - 
Since ks is very small (ks = 0.4), they are decreasing with v max . The characteristic non- 
monotonous behaviour for increasing k D (see is reproduced for all t> max values. Since 
in this parameter regime the particle density is a restricting factor for the evacuation 
time due to the clogging occuring at the exit, for t> max > 2 one finds some kind of 
convergence of the time values. 

4. Finer discretisation of space 

Increasing the maximal velocity of pedestrians is not the only modification that 
influences the realism of the model. One of the boldest assumptions in CA models for 
pedestrian motion is the representation of space as a grid of cells and the corresponding 
restriction of the spatial resolution to the cell size, in this case 40 cm. However, it can 
be argued that this simplification is justified by the reaction time (or decision time) 
which can be easily identified in the model. Nevertheless, an - preferably qualitative 
- investigation of the discretization effects would shed some light on the issue. The 
space discretisation used so far has been chosen in the simplest possible way such that a 
pedestrian occupies only one cell. This leads to a length scale of a = 40 cm (see Sec. EJ). 
However, for some applications this might not be adequate and a finer discretisation is 
necessary. In the following we will discuss the necessity of a smaller length scale, ways 
of its implementation and the consequences for the behaviour of the model. 

4-1. Motivation and consequences 

The identification of the size of one lattice site with the typical space occupied by 
a pedestrian in a dense crowd is most natural for the construction of a CA model for 
pedestrian dynamics. For square cells this leads to a cell size ofa 2 = l/p^ x p) and a = 40 
cm where pmSc P ' ) ~ 6.25 persons/m 2 is the maximal empirically observed density [14J. 
All units of space are then a multiple of the length scale a = 40 cm. In contrast to the 
introduction of higher velocities discussed in the previous section which was motivated 
by the unrealistic shape of the fundamental diagram, there are no simulation results 
that indicate the necessity of a better spatial resolution. Nevertheless there are good 
reasons to introduce a finer discretisation of space. 

(i) If the model is used not only for basic investigations of pedestrian dynamics, 
but also for the assessment of evacuation scenarios in complex structures (like 
passenger vessels or football stadiums), an accurate representation of the geometries 
is desirable. A finer discretisation corresponds to a more accurate representation 
of geometrical structures in a natural way. This can be seen from Fig. El 
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1 m 



Figure 9. Representation of a corridor of width 1 m in the 20-cm-model. Such a 
corridor can not be reproduced exactly in the 40-cm-model. The four cells occupied 
by each particle are shaded gray. 



(ii) As we have seen in Sec. El some situations require the use of higher velocities 
j) mM > 1. If one wants to keep the time scale unchanged this has to be compensated 
by introducing a smaller length scale. 

(iii) In principle it would be interesting to consider the continuum limit a — > in order 
to make contact with models that are based on a continuous representation of space, 
e.g. the social force fIE\ l2*§] and similar models [30J or hydrodynamic approaches 

jsnisaisaiiniiiSH!. 

The straightforward approach to investigate the discretization effects is to vary the 
cell size. If one chooses a cell size of a = 20 cm, the area A = 1/ p^f 1 occupied by a 
pedestrian is 4 cells (see Fig. EJ) in order to reproduce the maximal density pmax • The 
natural generalisation for the allowed motion of the original model as described in Sec. 121 
(in the following denoted as 40-cm-model) is the motion of the center of a pedestrian 
by one cell (i.e. a distance a, see Fig. 110)). This is only possible if the two neighbouring 
cells are unoccupied. As can be seen from the relation (JJJ) (with t> max = 1) for the three 
fundamental quantities the time step At in the 20-cm-model is half as long as in the 
40-cm-model^f. For the empirical velocity of t> emp = 1.3 m/s and a = 20 cm the time 
scale is At = 0.15 sec. For an arbitrary spatial discretisation a one obtains the time 
scale At = a/v^t . 




20 cm 



Figure 10. Occupation of 2 x 2 lattice sites of cell length a = 20 cm and possible 
transitions in a time step. 

If An extension of the 20-cm-model to a model with walking speeds u max > 1 is possible in the same 
way as discussed in the previous section, but will not be considered in the following. 



Discretisation effects and the influence of walking speed in cellular automata models for pedestrian dynamics. 



Finer discretisations of the underlying geometric structure for fixed velocity t> max 
therefore lead to smaller length and time scales of the model. So the limit a — > results 
in a space and time continuous model. This limit allows a comparison between discrete 
and continuous models, but is not considered here. 

The occupation of four lattice sites by one particle has further consequences in the 
20-cm-model, especially for the floor field model described in Sec. HJ 

(i) Every particle leaves two of the four cells it occupied previously (Fig. ITUj). In both 
these cells the dynamic floor field has to be increased. 

(ii) Since a particle occupies always four lattice sites, the transition probabilities are 
calculated by averaging the specific field values of static and dynamic floor fields S 
and D. 

In any CA, not only the floor field model, an important new effect related to conflicts 
occurs when the cell size is reduced and particles occupy more than just one cell. In the 
40-cm- model, conflicts can be described as local interactions (Fig. ITTk). because only a 
small number of particles can access the same lattice site which is determined by the 
coordination number of the lattice. On the other hand, in the 20-cm-model conflicts 
are not necessarily local anymore. Here, conflicts are not restricted to single cells and 
can spread over a wider area of the lattice. This is illustrated for a typical situation 
in Fig. fTTh. Assume that the permutation a (see Sec. IrS.l.lj) is such that the particles 
move from "left to right" , i.e. all the particles in the top row are allowed to move and 
the particles in the bottom row are not. Switching <t _1 (1) and <j~ 1 (2) will change the 
movement of all particles. This is not the case for particles that occupy only one cell, 
i.e. a = 40 cm. 

a) b) 




Figure 11. a) Typical local conflict in the 40-cm-model. b) In the 20-cm-model in 
principle all particles of a system can be part of a conflict. Hatched cells contribute to 
the non-local conflict. 

In the floor field model conflicts allow for the introduction of friction parameter /i 
[HI EH] (see Sec. I2J) • Due to the possible non-local nature of conflicts in the 20-cm-model 
this concept can not straightforwardly be generalized. New concepts for the description 
of conflicts are necessary, which is currently under investigation and will be discussed 
in a future publication. 
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In the next section simulation results of the 20-cm-model will be compared to results 
of the conventional model. 



4-2. Simulation results 

4-2.1. Evacuation times First, we will turn to the influence of the discretisation on 
evacuation times. We measured the averaged times of the evacuation of a large single 
room with an area of about (24.8 m) 2 , corresponding to 62 x 62 cells for the a = 40-cm- 
model or 124 x 124 cells for the a = 20-cm-model. The particles leave the room via a 
door of width 80 cm (two or four lattice sites, respectively). Fig. IT2l shows the impact 
of the variation of the parameters ks and krj in the a = 20-cm and a = 40-cm-model + . 
The initial particle density is p = 0.3 corresponding to 1153 particles. The qualitative 



a) 



b) 



k D =0.0, a=40 cm 
k n =0.0, a=20 cm 



2500 



k s =0.4, a=40 cm 
k a =0.4, a=20 cm 



Figure 12. Averaged evacuation times for a room of size 62 • 40 cm x 62 • 40 cm and 
an exit of width 80 cm: a) kr> = 0.0 fixed and variation of ks; b) ks = 0.4 fixed and 
variation of kn (5 = 0.3 and a = 0.0). 

behaviour of the model when varying the parameters (see Sec. Eland jjj) is not influenced 
by the finer discretisation. Because of slightly different transition probabilities due to 
the averaging of the field values of four sites for each walking direction in the a = 20- 
cm-model, small deviations only exist for coupling strengths ks < 3 (see Fig. 112b). This 
also leads to different time values for the parameter choice ko = 0.0 and ks = 0.4 in 
Fig. 112b). It is remarkable that the effect of the non-monotonic parameter dependence 
[7J of the curve is most pronounced in the a = 20-cm-model. 

The most important difference between the 20 cm- and 40 cm-model is the 
possibility of a deadlock in front of a small exit. Consider the flow through a bottleneck 
of the dimension of the particle width (i.e. bottleneck and particle have the width of 
two cells or 40 cm). If two particles reach the bottleneck at the same time step, their 

+ The times are measured in update steps. Therefore, for a comparison of the two models, the 
simulation times of the a = 20-cm-model have to be divided by two. 
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mutual blockage (see Fig. IT^j) may lead to vanishing flows. Similar observations have 
been reported in experimental investigations (see e.g. [SHIEZI)- This phenomenon can 

Exit 
A 




Figure 13. Blockage of bottlenecks or exits of width 40 cm in the 20-cm-model caused 
by two particles. Both particles try to access the four lattice sites in front of the exit. 

occur in pushing pedestrian crowds at bottleneck-shaped exits of small width. It has 
also been observed for some continuum pedestrian models with extended particles [2Zj 
and similarly in granular flow through a funnel . Therefore it is remarkable that 

it does not appear in the original 40-cm-model. 

On the other hand, the movement itself is more effective for the 20 cm-model 
compared to the 40 cm-model. For very high values of ks, the overall evacuation time is 
determined by the outflow through the door. In the case considered here (w = 80 cm = 4 
cells if a = 20 cm), the blockage does basically not occur. For very low values of 
ks (cf. Fig. I12j) . the particles perform a random walk and therefore the discretization 
(a = 20 cm vs. a = 40 cm) does not play a prominent role. This interpretation is also 
supported by the results shown in Fig. ITH1 where the smallest possible exit is investigated. 
In this case, the flow increases with the density for small coupling strengths ks and 
increases with large coupling strengths. However, this does not really explain the subtle 
differences in the area 0.5 < ks < 2, which will be investigated in more detail at a later 
point. 

4-2.2. Flow measurements The discretisation has a significant influence on the flow. 
Fig. El shows the fundamental diagrams in the ordered regime for a corridor of width 
2 m (i.e. width five lattice sites in the a = 40-cm-model and ten lattice sites in the 
a = 20-cm-model). The flows of both models are identical up to p ~ |. The density 
region p < \ characterises in both models the free flow regime. Nearly all particles can 
move to their right next neighbour site in each time step, the flow increases linearly 
with the density. For a density of p ~ ~, the probability of occupied next neighbour 
sites is highly increasing in the 40-cm-model, so the average velocity of the particles is 
decreasing to a value (v x ) < 1. Therefore, the flow decreases with increasing density 
and is additionally decreased by lane changing events (see Fig. HJ). 

The behaviour is different in the 20-cm-model: The space a particle needs for 
unimpeded movement in one time step corresponds to only half the length of its own 
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Figure 14. Fundamental diagrams in the ordered regime (ks = 10.0, kn = 0.0) of 
both model variants. Width of corridor: 5 cells in the 40-cm-model and 10 cells in the 
20-cm-model. 

size. Therefore, the average particle velocity is (v x ) « 1 up to a density of p ~ |, which 
leads to a linear increase of the flow up to this density. For p > | one finds a sharp 
decrease of the flow. 

Thus the finer discretisation results in a shift of the maximum flow to regions of 
higher density. This effect will be present for even finer discretisations. As in Sec. El 
experimental data point at a realistic flow maximum for densities p < |. Therefore the 
20-cm-model (with v max = 1) can not reproduce these results qualitatively. 

Fig. ITol shows the flow J through a bottleneck of width 40 cm versus the parameters 
ks and p. Only for very small ks {ks = 0.1) a small particle flow is sustained for higher 
densities. For increasing ks the flow breaks down at very small densities (for ks = 10.0 
the flow is zero for all densities). The system freezes for growing coupling strength ks to 
the static floor field. This is somewhat similar to the paradoxical f aster-is- slower- effect 
[2 |27] , as an increasing ks implies higher individual walking speeds. The clogging effect 
observed in our model is due to the local blockage shown in Fig. [T31 

5. Conclusion 

We have systematically investigated the influence of the interaction range v max and the 
spatial discretisation a on the behaviour of two-dimensional CA models for pedestrian 
motion. Using simulations of a corridor and evacuation from a large room we have shown 
that both parameters can have a significant influence on the properties of the model. 
Although we have used a specific model, the floor field CA, for our simulations due to 
its ability to reproduce the observed collective phenomena correctly we believe that the 
effects discussed here are generic for any discrete approach to pedestrian dynamics. Of 
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Figure 15. Fundamental diagram of a corridor (93 x 33 cells of size 40 cm x 40 cm) 
with a bottleneck of width 40 cm. The couplings to the static and dynamic floor fields 
are k$ and ku = 0.0. 

course this requires further investigations in the future. 

First we have considered the extension of the interaction range which we have 
implemented through a larger walking speed v max . We have seen that there are different 
possibilities to achieve this. This modification leads to a considerable improvement of 
the agreement of fundamental diagrams for the motion along a corridor with empirical 
findings. In the original case t> max = 1 the fundamental diagram of the model was 
almost symmetric with respect to the density p = 0.5. The reason is that for a such 
a uni-directional motion where swaying is strongly surpressed many models effectively 
reduce to a number of independent ASEPs* and hence the observed flow-density relation 
reflects its particle-hole symmetry. As in models of highway traffic J7j the maximum 
in the homogeneous case is shifted to lower densities for larger t> max = 2 and thus leads 
to a more realistic approach. 

As a second scenario we have studied the evacuation from a large room. Here it 
turned out that in most cases the results are not affected strongly by the introduction 
of a larger velocity v max . 

Summarizing, these results indicate that the choice w max = 1 might not in all cases 
be justified. As indicated above the fact that for the corridor one has to deal with a 
basically one-dimensional motion is important. Here it is remarkable that also for a 
continuum model like the social force model this problem requires a rather different 
choice of interaction parameters [37] compared e.g. to evacuation problems, especially 
concerning the repulsive interaction forces. This shows that both model classes have 
to be generalized in order to achieve a proper description of all possible situations in a 

* One could expect that in the case of a mixture of species with different v m&x the decoupling of the 
lanes will be less pronounced. 
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unified way. 

The case of the second non-trivial fundamental parameter, namely the cell size, is 
even more intricate. A finer spatial discretisation does not only provide a connection 
to a different class of microscopic models based on continuous representation of space. 
In many applications to real life problems the geometries under consideration can not 
be described by lengths that are multiples of the length scale a = 40 cm of the original 
model. We have argued that a reduction of this length scale, and thus the cell size, 
must be accompanied by an increase of the size of the particles. In order to reproduce 
the empirically observed maximal density they have to occupy now more than one cell, 
e.g. four in the case of a = 20 cm. If the particles should have the same symmetry as 
the underlying lattice this also restricts the allowed length scales. E.g. for the square 
lattice the most natural extensions of the particles are 2 n x 2 n cells corresponding to a 
length scale a n = ^ cm. 

By studying the same two scenarios as for the influence of t> max we found two main 
consequences of using a finer discretisation. On the one hand, the flow is increased since 
the particles can now move ahead if there is space available less then their own size. 
On the other hand, particles can now block each other more easily. This might lead 
to non-local conflicts where groups of pedestrians mutually block each others motion 
(Fig. n"Tj) . Furthermore, at exits blocking can result in a complete breakdown of the 
outflow. This is not related to the discreteness of the model and can also be observed 
in continuum models with particles of finite size. 

It would be interesting to perform the continuum limit of our model to relate it 
to other continuum approaches. A naive continuum version of the floor field model is 
currently under investigation and might serve as a reference model. Our present study 
has shown that finer discretisation or larger velocities in combination with a discrete 
time (parallel) dynamics leads to the possibility of complicated conflict situations. One 
could try to solve this problem by using a continuous time dynamics, but this might 
drastically increase the computational complexity of the model in a different way. 

Even though our investigations are based on a square grid (one could also think 
of a hexagonal lattice, and indeed there are models based on this [4*0*1 141 j ) the results 
obtained are generally valid. As a next step currently the combined effect of larger v max 
and smaller cells is investigated. Then one can expect to encounter a combination of 
many effects observed here. However, e.g. in the case of flow along a corridor these 
modifications have opposite quantitative effects. Whereas an increase of t> max tends to 
shift the location of maximal flow to smaller densities the reduction of the length scale 
has the opposite effect. This work is in progress and will be reported elsewhere. 

Summarizing we can say that discrete models like cellular automata yield quite 
realistic results. E.g. the floor field model is able to reproduce collective effects (such 
as lane formation in counterflow etc.) observed empirically. However, for applications 
where also a reliable quantitative prediction should be made one has to take some care 
in the choice of the basic parameters as cell size and interaction range. 
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